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ABSTRACT 


This report presents the theory, computer program and numerical 
results for an axial slot antenna on a circular cylinder. 
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I. INTRODUCTION 


We consider an axial slot antenna on a perfectly conducting circular 
cylinder. The cylinder is partially coated with a dielectric layer, and 
the antenna radiates through this flush -mounted window. The motivation 
for this study is to determine the effects of a high-temperature dielectric 
layer on the performance of antennas mounted on a space shuttle. 

For an axial slot antenna on a circular cylinder completely coated 
with a dielectric layer, the admittance and patterns have been investigated 
by Knop[l], Fan,te[2], and Croswell, Westrick and Knop[3]. Our analysis has 
some similarity to that of Billingsley and S1nclair[4] for scattering by 
circular-sector cylinders. 

The following sections define the problem and present the theory, 
computer programs and some numerical results. 


II. THEORY 

Consider an axial slot antenna on a perfectly conducting circular 
cylinder as illustrated in Fig. 1. The inner aperture has radius “a" and 



Fig. 1. An axial -slot antenna radiates through a flush-mounted 
dielectric window in a conducting circular cylinder. 
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half-angle The outer aperture has radius b and half-angle The 
exterior medium is free space. The inner slot radiates through a flush- 
mounted homogeneous dielectric window with permittivity e] , permeability 
s inner radius a, outer radius b and half-angle The metallic flange 
prevents the dielectric window from falling out. This cylindrical structure 
has infinite length, and its axis coincides with the z axis. We consider 
a time-harmonic excitation with the time dependence understood, and 
the fields have no z dependence. This report considers the TE polarization 
in which the non-zero field components are Ep, E^ and H 2 . Given an even 
field distribution over the inner aperture, the objective is to deter- 
mine the aperture admittance, gain and far-field pattern of this antenna. 

Our solution employs cylindrical -mode expansions and Galerkin's method. 


The field in region I (the dielectric window) is 


(1) 

E>. 

P 

k] p 

V [C|^ J^(p) 

+ \ sinv(|> 

(2) 


I 

‘ k 

j;{p) + 

d|^ N^(p)] cosvi^ 

(3) 


I [C|^ 

k 

Jy(p) + d,^ 1 

Ny(p)] cosv<() 

(4) 

h * 

e-| 

- 


(5) 

"1 ■ 




(6) 

V = 





where the integer k runs from zero to infinity and (p,<J>,z) are the 
cylindrical coordinates. (In this report the symbols Jv(p) and N^(p) 
denote the Bessel and Neumann functions with order v and argument k]p.) 
This field satisfies the source-free version of Maxwell's equations in 
region I. From Eqs. (1) and (6), tangential E vanishes at the perfectly 
conducting surfaces at ^ = ± (J> 5 . The expansion constants C|^ and d|^ are 
to be determined from the boundary conditions. 

The voltage across the inner aperture is 
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( 7 ) 


V = 2a E 
'o 


f: 


where E denotes £ 4 ( 3 , i|>). (We assume the aperture field E is a specified 

even function of $.) The external admittance of the inner aperture is 

( 8 ) 


Y = 


- 2a 
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] E* H^(a,i))d<|> . 


The boundary condition at p = a is 

E for 0 < (j) < (j). 


(9) 




0 for <|ij^ < i}* < <J>jj 


From Eqs. (2) and (9) with Fourier analysis, 
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= E cosv<ti d<|> 
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where e^ = 1 and = 2 for k =1, 2, 3*** 

From Eqs. (3) and (8), the external adnittance (per unit length of cylinder) 
of the inner aperture is 

(12) Y = ^ HC|( J^(a) + d|^ N^(a)] G* . 

k 

The field in region II (the exterior free-space region) is 
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where the integer i runs from zero to infinity. (In this report the symbol 
H-j(p) denotes the Hankel function with order i and argument I<qP and the 
superscript (2) is understood. The argument k-]p will not be encountered 
with the Hankel function.) This field satisfies the radiation conditions 
and the source-free version of Maxwell's equations. 

To complete the solution, it remains only to enforce the boundary 
conditions at p = b. The rigorous solution involves an infinite system 
of simultaneous linear equations. We desire an accurate approximation 
involving a finite system of simultaneous linear equations. To develop 
a solution of this type, we expand the field in the outer aperture 
(at p = b) as follows: 


(18) E. = I cos(mr(|)/(|)^) for 0 < (j> < 


where n runs from zero to N. If the constants b^ were known, the remaining 
constants (ai, c|< and d^) could be determined. In this sense the bn are 
independent unknowns, and the others are dependent. When the simultaneous 
linear equations are written as a matrix equation, the square matrix will 
be symmetric if the bn are chosen as the independent quantities. 

I 

From Eq. (2) and the boundary condition on at p = b, 


(19) 


jni I [c,^ d^(b) + d|^ N^(b)] cosv<j. = 


for 0 < <t> < 

(j) c 

0 for (j)^ < 


where E. is defined by Eq. 
cosv<t> and integrating over 


(18). Multiplying both sides of Eq. (19) by 
the range 0 < (|) < <|)|^ yields 


(20) Jn, [Ck * “k J “’n 
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(21) 
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From Eqs. (10) and (20), 
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From Eq. (14) and the boundary condition on at p = b, 


(25) 


3% I H! (b) cos i<|> = 


E^ for 0 < 4> < <|>^ 
0 for < <fi < n 


where E^ is defined by Eq. (18). Multiplying both sides of Eq. (25) by 
cos i+ and integrating over the range 0 < ()> < tt yields: 
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Equations (22)-(26) show explicitly that a knowledge of the constants 
bfi is sufficient to determine all the other constants. 

At this point we have used the boundary conditions oh E^ to relate 
a-j, C|(; and d^ to b^. The next step is to use the boundary condition on 
to generate a system of simultaneous linear equations for the constants 
bp. From Eqs. (3) and (15) and continuity of tangential H across the 
outer aperture (at p = b): 
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( 28 ) 


J (b) cos i«|) = J [cj^ J^(b) + d^^ N^(b)] cosv<f. 

1 k 

where <f> ranges from zero to <l)c. In Eq. (28), multiplying both sides by 
cos(niiT(j)/(|)(,) and integrating from <{> = 0 to yields 


(29) 
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In matching across the aperture, we selected the same weighting 
function cos{imr(j)/(J»(.) also used as a basis function in Eq. (18). This is 
the distinctive feature of Galerkin's method. If Eqs. (22), (23) and (26) 
are used to eliminate a^ , and d^, Eq. (29) yields: 
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Equation (30) is recognized as a system of simultaneous linear equations. 

In the summation, n runs from zero to N. Equation (30) can also be 
written as a matrix equation. The symmetry of the square matrix is 
obvious in Eq. (32). 

The matrix equation is solved with a digital computer to obtain 
nunrerical values for bp. Then Eqs. (22), ^3) and (26) are employed to 
determine ci., dj^ and aj. The aperture admittance is obtained from Eq. (12), 
The far-field pattern is obtained from Eq. (15) as follows: 


(34) 




cos iifi . 
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The power gain is calculated as follows: 


(35) 


Gp(4») 


2irpn iH (p,<f)!^ 


where the aperture voltage V is given by Eq. (7) and the aperture con- 
ductance G is the real part of Y in Eq. (12). 

III. NUMERICAL RESULTS 


Figures 2-6 illustrate the far-field patterns of an axial slot 
antenna radiating through a lossless window with dielectric constant of 
1, 1.2, 2, 3 and 4. Figures 2 and 3 compare the calculated patterns with 
experimental measurements performed at NASA Langley. Measurements are 
not available for the other cases. In this sequence of figures all 
parameters of the slot, window and cylinder are fixed except the dielectric 
constant. The electric field distribution is uniform across the inner 
aperture. 


All the patterns are reasonably smooth except in Fig. 5 where the 
pattern breaks up into many lobes with deep nulls. With a dielectric 
constant of 3, this anomalous type of pattern is observed when the aperture 
half-angle is <|>|j = 13.8, 14.8, 15.8, 16.8°, etc. When differs from one 
of these critical angles by more than 0.1 degrees, the pattern becomes 
smooth again. At each critical angle, the aperture width is an integral 
number of wavelengths for the lowest-order surface wave. This surface- 
wave resonance phenomenon is less pronounced with a dielectric constant 
of 1.2 but is observed when <t>k = 14.6°. The effect may be reduced with 
a lossy dielectric window or by reducing the reflection coefficient at 
the edges of the aperture. 

In these figures, the calculations are based on a two-dimensional 
model with an infinitely long axial slot. A case of greater interest is 
a half-wave axial slot in a long cylinder. The effects of surface-wave 
resonance will be reduced with a slot of finite length. 

In generating the data for Fig. 5, the execution time was 50 seconds 
on a Datacraft 6024/3 computer. The solution involved a system of 20 
simultaneous linear equations (N = 19 in Eqs. (18) and (30)). The 
infinite series with index i (in EqS. (32) and (34)) were truncated after 
148 terms, and the series with index k (in Eqs. (12), (31) and (32)) were 
truncated after 20 terms. The calculated aperture admittance was 0.372 + 
j 0.122 mhos/wavelength. Identical results were obtained with N = 18 and 
19, but N = 15 proved inadequate. 
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Fig. 2. Far-field pattern with = 1 . 
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CALCULATED 


— MEASURED 


Window permittivity: 


= 1.2e 


Inner radius: 
Outer radius: 


18.7325x^ 

19.05x„ 


Angular extent of inner aperture: ~ 1.08' 


Angular extent of dielectric cavity; 
Angular extent of outer aperture: 


2 ^^ = 29.6“ 
2<)>^ = 29.6° 
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Fig. 3- Far-field pattern with Cy, = 1.2, 
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Fig. 4. Far-rfield pattern with = 2 . 
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Fig. 6. Far-field pattern with Cp = 4. 
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IV. SUMMARY AND CONCLUSIONS 

This report develops the theoretical formulation for a TE axial 
slot antenna radiating through a dielectric window in a circular cylinder. 
Numerical results are presented for the far-field patterns, and it is 
noted that the calculations show excellent agreement with experimental 
measurements. The computer program is presented in the Appendices. 

The solution is based on Galerkin's method. Simultaneous linear 
equations are generated in which the unknown quantities are the coef- 
ficients in a Fourier-series expansion for the electric field in the 
outer aperture. The formulation is rapidly convergent, and the computer 
program is quite efficient. 
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APPENDIX I 

THE MAIN COMPUTER PROGRAM 


The MAIN computer program is listed in Fig. 7. In this program 
is uniform across the inner aperture and the dielectric window is 
lossless. Following the format statements, the dimensions are Indicated 
for the subscripted quantities as follows: 

IDC dimension of B and V 

IDH dimension of A, BHR, BB, YY, BP, YP, SNC and SGA 
IDJ dimension of C, D, SGC, RBES, BEN, AJJ, etc 
IDZ dimension of Z. 

The input data are prograrrmed at statement 20 with the following definitions 

AL inner radius a/x 

BL outer radius b/x 

ER dielectric constant e-j/eo 

DPH angular increment for far-field pattern calculations 
PHA <j>g, in degrees 

PHB (j)jj in degrees 

PHC (jij. in degrees 

where X denotes the wavelength in free space. 

At statement 30, subroutine BESSI is called for the Bessel and 
Neumann functions and their derivatives. This subroutine also determines 
the number of terms (denoted by KK) to be employed in the summations on k 
in Eqs, (12), (31) and (32). The last call to BESSI determines the 
number of terms (denoted by II) to be employed in the summations on i in 
Eqs. (32) and (34). If KK is equal to IDJ, the dimension IDJ should be 
increased. If II is equal to IDH, the dimension IDH should be increased. 

II should exceed IMIN, and KK should exceed KMIN. NN denotes the number 
Of simultaneous linear equations and the number of terms to be employed 
in the summations on n in Eqs. (22), (23) and (26). NN should exceed 
NMIN. 

For a uniform aperture distribution with V = 1 volt, Eqs. (7) and 
(11) yield 

sin(v^_) 

ol 

In the computer program, SGA(K) denotes 2aG|(. Subroutine GLJ calculates 
Ekn/‘('c where F^n is defined by Eq. (21). Subroutine GNJ calculates 
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G-in/<j»c where G-jn is defined by Eq. (27), Some of the symbols used in the 
program are defined as follows: 


A 

B 

C 

D 

V 

AK 

BK 

AKl 

BKl 

BEN 

BHR 

ETA 

ETAl 

GNU 

RBES 

S6C 

SNC 

SJN 

SJL 

Yll 

Z(L) 




“k 

k-|d 

k^b 


denominator in Eq. (33) 
H. (b)/H](b) 


'0 

^1 

V 


(sinvi>^)/(v(j)^) 

(sin i<(>^)/(i«j>^.) 
first summation in Eq. (32) 
second sunination in Eq. (32) 
aperture admittance Y 

mn 


In statement 130, subroutine SQROT is called to solve the system 
of simultaneous linear equations. Then the expansion coefficients a-j , 
cu and \ are calculated from Eqs. (22),. (23) and (26). The aperture 
admittance is calculated at statement 360 with Eq. (12). Finally, the 
gain is calculated with Eqs. (34) and (35). 
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TE AXIAL SLOT IN PERFECTLY CONDUCTING CIRCULAR CYLINDER. 0001 

SLOT RADIATES THROUGH DlrLECTRIC WINDOW. 0002 

PROGRAM BY JACK H RICHMOND, OHIO STATE UNIVERSITY. 0003 

COMPLEX vr 30I,B( 30),2( A65),C( 150J,0( ISO ) , A t 400 » , HHK ( 400 ) 0004 

complex CO,HZ,SUMM,SUL,SJN,Yll 0005 

DIMENSION B»(400),YY(400I,BP( 400 I , YP ( 400 I , SNC ( 400 ) , SGA ( 400 ) 0006 

DIMENSION SGCI 150) ,R»ES( 1501, BENI 150) .0007 

DIMENSION 4JJ( 150), AYYI 15()),AJP( 150),AYP< 150) 0008 

DIMENSION BJlt 150 ) ,BY1 ( 150 ),RJPl( 150) ,BYP1( 150) 0009 

equivalence IB, V) , (BB, SNC ), ( YY,SGA), (BJl ,SGC) , I BYURBES) 0010 

DATA ETA, PI, TP/376. 727, 3. 14159, 6. 28318/ 0011 

2 FORMAT! IX, BF15.6) 0012 

4 FORMAT! IX, 121 10) 0013 

5 FORMAT! IHO) 0014 

IDC=30 0015 

10H=400 0016 

IDJ=150 0017 

I0Z=465 0018 

20 AL=18.7325 0019 

BL=19.05 0020 

ER=3. 0021 

0PH=2. 0022 

PHA=0.54 0023 

PHB=14.8 0024 

PHC=PHB 0025 

IF! PHA.GT. PHB)PHA=PHB 0026 

IF! PHC.GT,PHB)PHC=PHB 0027 

N=.5+(S0RT( 1.+8.4IDZ )-l. )/2. 0028 

IFIN.LT. IDC ) IDC=N 0029 

NN=IOC 0030 

SOR=SORT(ER) 0031 

£TA1=ETA/S0R 0032 

TL=BL-AL 0033 

WRITEI6,2)AL,BL,TL,ER,PHA,PHB,PHC 0034 

WR1T£!6,5) 0035 

AK=TP4AL 0036 

BK=TP*BL 0037 

AK1=AK*S0R 0038 

BK1=BK«S0R 0039 

PHAR=.bl74533*PHA 0040 

PHBR=,0174533*PHB 0041 

PHCR=. 0174533=»PHC 0042 

GNU=PI/PHBR 0043 

30 CALL BESSI! AK1.GNU,AJJ,AYY,AJP,AYP, IDJ, 1DH,KK ,BB,YY) 0044 

CALL BESSI I BK1,GNU,BJ1,BY1,BJP1,BYP1, IDJ, I0H,LL ,BB,YY) 0045 

CALL BESSI ! BK, 1. ,68,YY, BP, YP, IDH, IDH, I I ,BB,YY) 0046 

1F(LL.LT.KK)KK=LL 0047 

KMIN=BK4PHB/180. 0048 

IMIN=BK 0049 

NMIN=BK4PHC/180. 0050 

NMAX=.5+KK«PHC/PHB 0051 

IFINNiGT. IOC)NN=IOC 0052 

1F(NN,GT,NMAX )NN=NMAX 0053 

IFINN.LT. 1 )NN=1 0054 

WR ITE 16,4 ) IMIN, 1 1 ,KMIN,KK,NMIN,NN 0055 

WRITE(6,5) 0056 

SUMN=(.0,.0) 0057 

BHR( 1 )=CMPLX| BB! 1 ),-YY( 1) I/CMPLXIBPI 1 ),-YP! 1 ) ) 0056 

00 60 1=2,11 0059 

BHRI I )=CMPLXI BB! I ),-YY! I ) l/CMPLX!BP! I ),-YP( I ) ) 0060 

N=I-1 0061 

SC=S1N(N4PHCR)/!N*PHCR) 0062 


Fig. 7. The MAIN computer program. 
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SNCU)*SC 0063 

60 SUMN = SUMN+BHR( n*SC*SC OO*'^ 

RUM1=AJP( I >-6Yl( 1 )-hJH n#AYPl n 0065 

REN1=AJP( n*BYPl( n-BJPl< 1 j^AYPU ) 0066 

SUML=.0 0067 

SUMP=.0 0068 

DO 70 K=2,KK 0069 

BEN(K»=AJP(K>«BYP1(KJ-BJP1(K)*AYP<K» 0070 

RBES(K) = ( AJP(K)*RYl(K»-BJUK)«AVP( K) )/BEN(KJ 0071 

GNU=(K-1. )*PI/PHBR 0072 

SC=S I N( GNU>»PHCR ) / ( GNU^PHCR » 0073 

SGC{K)=SC 0074 

SGA(K)=S!N(GNU*PHAR)/(GNU«PHAR» 0075 

SUMP=5UMP+SGA(K)*SC/BeN(K) 0076 

70 SUML=SUML + RBES< K)=i‘SC’i‘SC 0077 

Z( 1)=ETA1*PHBR*«.5hcBHR( 1)+SUMN)/(ETAX'PI 1-.5<'RUMI/RENI-SUML 0076 

2 ( n=PHCR>»Z< 1 ) 0079 

V ( 1 »=CMPLX( < l./RENl+2.*SUMP)/( AK^BKl ),0. ) 0080 

IFlNN.EO.nCO TO 130 0081 

00 120 N=l,NN 0082 

MA=2 0083 

IF(N.GT.21MA=N 0084 

DO 120 M=MA,NN 0085 

SJN=(.Or.O) 0086 

00 100 1=2,11 0087 

CALL GNJ«I,M,PHCR,SNC,GIM) 0088 

GlN=SNCnt 0089 

IF(N.GT, 1 )CALL 6NJ( 1,N,PHCR,SNC,GIN) 0090 

iqO SJN=SJN+BHR( I )>t=GIM*GIN 0091 

SJL=«0 0092 

CJR=.0 0093 

00 110 K = 2,ICK 0094 

CALL GLJ(K,M,PHBR,PHCR,SGC,FKM) 0095 

FKN=SGC(K) 0096 

IFlN.GT.nCALL GL J ( K ,N, PHBR, PHCR , SGC, FKN ) 0097 

IF(N,EO. 1 )CJR=CJR+FKM«SGA(K»/6EN(K) 0098 

110 SJL=SJL+RBeS(K»*FKM*FKN 0099 

L=(N-1 )<‘NN-tN*N-N)/2 + M 0100 

Z(L)=<ETA1*PH6R«SJN/<ETA*P1)-SJL)*PHCR 0101 

120 IFIN.EO. 1 )V(M)=CMPLX(2.*CJR/( AK*8K1 ),0. ) 0102 

130 CALL SOROT( Z,V,0, 1 ,NN) 0103 

A( !)=-( . 0 , 1 . )*PHCR*B( I )/<ETA«PI*CMPLX(BP(l l,-YP(im 0104 

DO 300 1=2,11 0105 

SC»SNCm 0106 

SUL*Bm»SC 0107 

IFlNN.EO. 1160 TO 300 0108 

DO 290 N=2,NN 0109 

CALL GNJ( I,N, PHCR,SNC,G1N1 0110 

290 SULsSUL+BlNlfGIN 0111 

300 Ad )=-2.*PHCR«( .0, 1. )*SUL/( ET4*P I^CHPLX < BP ( I),-YP( I ) >| 0112 

Cl 1 ) = ( .0, 1, 1*( 2.*AL»PHCR<'B( 1 )*AYPm-BYPl( I ) ) / 0113 

2<2,*AL#ETAl*PHBR>KReNl) 0114 

0m = ( .0, 1. )*(BJP1 ( 1 )-2,fAL*PHCR*Bm'*AJP( 1))/ 0115 

21 Z.’O'AL^ETAl^^PHBR^RENl ) 0116 

DO 340 K=2,KK 0117 

SA=SGA1K) 0118 

REN = AL*ETA1*PHBR“1‘BEN(K) 0119 

SUL=ei 1 K'SGCIKI 0120 

IFlNN.EO. UGO TO 330 0121 

DO 320 N=2,NN 0122 

CALL GLJ1K,N, PHBR, PHCR, SGC, FKNI 0123 

320 SUL = SUL-»-BlN)«FKN 0124 


Fig. 7. 
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330 C(K*=-'SA*BYPHK H-2.^AL^PHCR-AYP(K)^SUL 
C( K »=< .0, 1, )«C(K)/KEN 

OIK l-SA’^BJPnK)^2,*AL=i‘PHCR*AJP{ K)«SUL 
340 Of K) =1 *0, I. J^O(K)/REN 

Y11=C( 1 J^AJJI 1 )-»-0( i )«AYY( I > 
on 360 K=2,KK 

360 Yll^YU + ICIK )’^AJJIKKLHK)«AYY( K ) )*SGA| K) 
WRITE(6,2)Yll 
WRITEI6t3> 

GG=:RFALI Yll ) 

NPH=IB0./0PH + K5 
DO 400 L=1,NPH 
PH=( L'-1)=*=DPH 
PHR=,0174533*PH 
HZ=^I .0,.0) 

C0^(U,0.) 

DO 390 I=lt I I 
N-r-1 

HZ=HZ+CO-A| 1 )«COSIN*PHR) 

390 CQ-CO>;‘l ,0t !• ) 

HAB==CABSCHZ)/Pl 
GAIN^TP’5<ETA*HAB^-HAB/GG 
0B = 10*=4'AL0G10IGAIN) 

400 WRITE(6,2)PH,GAlNtOB 
CALL EXIT 
END 
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APPENDIX II 
SUBROUTINE BESSI 


The subroutines are listed in Figs. 8-11. BESSI is a drastically 
modified and streamlined version of a program developed by Nelson Ma of 
the Department of Engineering Mechanics at The Ohio State University. 

This program calculates Bessel and Neumann functions and their derivatives. 
The argument must be positive and real. The order is positive and real, 
and it may be integer or noninteger. For the gamma function, BESSI calls 
subroutine GAMMA from the IBM 360 scientific subroutine package. The 
input data are defined as follows: 

X argument, greater than zero 
ORD order, greater than zero 
IDL dimension of BJJ, BYY, BJP and BYP 
IDM dimension of BJ and BY 

BJ and BY are work arrays for internal use. If ORD is an integer, BJ 
and BY may have the same names in the calling program as BJJ and BYY to 
reduce storage requirements. This is illustrated in the third call to 
BESSI in Fig. 7, The output data are defined as follows: 

BJJ(I) J^(x) with I = 1, 2, 3, ••• N and v= (I - 1 )*0RD 
BYY(I) N^(x) 

BJP(I) jJ|(x) 

BYP(I) n;(x) 

N maximum value of I 


N will not exceed IDL. If IDL and IDM are sufficiently large, N will be 
determined by the condition that BJJ(N) is less than 10-B or BYY(N) is 
greater than 10°. Comparison with other subroutines indicates that the 
output of BESSI may be accurate even when x is as large as 2000. The 
upper limit on x is not known. 

One call to BESSI generates a series of Bessel and Neumann functions 
with different orders. For example, if ORD = 0.5 the functions will have 
orders 0, 1/2, 2/2, 3/2, 4/2, 5/2, etc. If ORD = 2, the functions will 
have orders 0, 2, 4, 6, etc. These are the orders required in boundary- 
value problems involving wedges and circular-sector cylinders. 

BESSI uses the recursion techniques of Reference [5]. For x greater 
than 10, the phase amplitude-method is employed[6]. 
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In line 39, the user should replace l.E-38 with the smallest number 
his computer can handle without underflow. To obtain a few more Bessel 
and Neumann functions in the series, one may replace l.E-6 with a smaller 
number in line 69, and replace 1.E6 with a larger number in line 154. 

To obtain the maximum available number of Bessel and Neumann 
functions in the series, the required dimensions may be estimated as 
follows when x is greater^lthan one: 

IDM = 1.2 X + 100 - 1500/(x + 20) 

IDL = IDM/ORD, 



SUBKOUTJNE BESS U X ,OROt B J J , BYV , B JP , HYP , lOL t ! OM, N , BJ , BY ) 1 

DIMENSION Rjr I )fBY( 1 I, BJJU ), BYY( 1 1 ),8YP{ 1 J • 2 

DATA A, PI/. 577215665, 3. 3 
DATA CO,C1,C2,C3,C4,C5,C6,C7,C8,C5,C10,C11,C12,C13,C14,C15,C16 4 

B, C17,C18,C)9,C20,C21,r,22,C23,C24,C25,C26,C27/ 5 

C. 25,. 15625, -.375,. 1171875,-1. 15626, 1 . 8 75 , . 952 1 48436 E- 1 , 6 

0-2.36671875, 14. 2265625, -19. 68 75, -. 809 326 1 72E-1 , -4. 10058593, 7 

E5B.2246094,-277.875,354.375, .4l6666667E-l,-.25,.0125,-.35, 8 

F. 5580357186-3, -.4241 07143, 3. 60267857, -5. 625 ,. 3036 1944E-2 , 9 

G-. 486111111, 10. 2864583, -58. ,78. 75/ 10 

J=0 ll 

IFIX.LE.O. IGD TO 1 12 

IFIDRO.LE.O. »G0 TO I 13 

GO TO 2 1*^ 

1 N=J-1 15 

RETURN 16 

2 EA=2./X 17 

INT=0RD+,5 18 

1N=1000.*(ORO-1NT 1 19 

TLOG=ALar,( X/2. ) 20 

PIH=2./PI 21 

T2=1./(X#X» 22 

P14=4,/PI 23 

GAMMIsPIHK) A+TLOG) 24 

KMAX=X+10.«( 2. ■i'X#*. 333333+1, ) 25 

S0PX=S0RT(.5*PI*X) 26 

TPXsa./IPI'i'X) 27 

10 J=J+l 28 

JM=J-1 29 

FNUP=JM«tiRU 30 

N=FNUP 31 

FNU=FNUP-N 32 

IF( IN.EO.OlFNUs.O 33 

IF(IN,E0.01N=1 34 

NP1=N+1 35 

IFINPI.GT. lOMlGO TO 1 36 

NM1=N-1 37 

KsKMAX 38 

IFIK.LT.NPl .AND. IN.NE.OIGD TO 1 39 

1«:K 40 

ejc=.o 41 

BJB=l,E-38 42 

EB=EA>M 1+FNU) 43 

35 BJA=6B*BJB-BJC 44 

IFd.UE.lDMlBjn )=BJA 45 

EB=EB-EA 46 

BJC=BJB 47 

BJB=6JA 48 

lsl-1 49 

IF( l.GE. 1 )G0 TO 35 50 

IFIK.GT, 10M)K=I0M 51 

M = IK-n/2 52 

IF(X.GE.10.»G0 TO 59 53 

PHI=FNU+2. 54 

H0=3 55 

ALF = PHI>KBJ( 3)+Bjn) 56 

on 39 1=2, M 57 

M0=M0+2 58 

FM2=2’CI 59 

FM1=I-1 60 

F1=I 61 

7EMP=( I FNU+FM2)*|FNU+FM1 ) ) / t F 1* I FNU+FM2-2 . 0 ) 1 *PH1 62 


Fig. 8. Subroutine BESS I. 
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PHI=TeMP 63 

39 ALF=PHf«RJ( M0I4ALF 64 

GAMM=GA«MAtFNU + l. 1 65 

ALF=FA*«iFNU>»GAMM>!<ALF 66 

41 AJl=l. 67 

AJ2=1. 68 

JAN=0 69 

RALF=1./ALF 70 

00 43 1 = 1, K ““ 71 

IFUAN.FO.UGO TO 43 72 

IFr AJ1 .lt. l.F-6 .AMO. AJ2, L T. 1 . F-6 ) JAI\t=l 73 

BJ( U=ejn 14RALF 74 

AJ1=AJ2 75 

AJ2=ABS(8J( n ) 76 

1MX=! 77 

43 CONTINUE 78 

K=IMX 79 

M=(K-n/2 80 

IFIIN.NE.O .AND. IMX.LT.NPUGO TO 1- 81 

60 TO 100 82 

59 KOl»NT = l 83 

GNU=FNU 84 

61 ALl=GNU**2-.25 85 

A2 =C0*AL1 86 

A4 = (C1»AL1+C2)>«<AU 87 

A6=nC34ALl+C4)*ALl+C5 MALI 88 

A8=n (C6=!=AL1+C7)4AL1+C8)«AL1+C9MAL1 89 

A10=( M (C10«AL 1+CIHWAL1+C12)*AL1+C13)*AL1+C14)>!-AL1 90 

B=n ( A10*T2 + A8M=T2+A6 MT2+A4 MT2 + A2 91 

BNU=B=>T2+1.0 92 

ANU=BNU/SOPX 93 

A2=.5>»AL1 94 

A4=tC15*ALl+C16 MALI 95 

A6=HC17>*AL1+C18»#AL1+.75J>!‘AL1 96 

A8 = n (C19*ALl+C20)=i‘ALl+C21 MALI+C22)»AL1 97 

A10=( M (C23*AL 1+C24 )» AL 1+C25 M AL 1 +C26)>i!AL I +C27 MAL 1 98 

B=( (( A10>»T2+A8»«T2+A6)*T2 + A4)<-T2+A2 99 

TPHMB»*'T2+1.0 100 

PHl = TPHMX-(GNU + ,5)/P|H 101 

F1 = ANU«C0S( PH! ) 102 

Yl=ANUfSIN< PHI! 103 

IF(K0UNT.6T. 1 )G0 TO 65 104 

FSAVE=F1 105 

BY!l)=Vl 106 

6NU=FNU+1.0 107 

KDUMT=2 108 

GO TO 61 109 

65 F2=F1 110 

BY{2)=Y1 111 

F1=FSAVE 112 

ALF=BJ(2>/F2 113 

IF! ABSJFl ).GT.ABS(F2) )ALF = BJ(1 )/Fl 114 

GO TO 41 115 

100 1F<X.GE. 10. )GC) TO 150 116 

ARG=FN()>t'P I 117 

GARG=GAMM*#2 118 

IF(FNU.EO.O. (GO TO 116 119 

TERM=( 1. /PI (^EA^^tZ.fFNU) 120 

6AM1=C0S(ARG)/SIN( ARG)-TeRM«(GARG/FNU! 121 

GAM2 = 2.0«TERM4GARG>1‘( FNU+2.0)/( l.O-FNU) 122 

GO TO 117 123 

116 GAMMGAMMl 124 


Fig. 8. 
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GAM2=PI4 125 

117 RY(2)=-( K/FI +2.«FNU)*GARG+(GAMI^GAM?/2p 126 

YNU-GAM1*BJ( 1 ) 127 


TXNU = 3.0={‘FNU/X 
AB=^ABS(BJf 1 ) ) -Op 00000^ 

MPl=M+i 
12 = 1 

on 121 I'2 tMP1 
12=12+2 
FI^I 
FIM=I-l 
FI2=2«I 

DEN0M = F I#( F 1-FNU ) ( FNl» + F I 2-2 p 0 ) 

' GAM3=(FNi; + F I 2 ) ={' { 2 . 0=^ FNU + F I M ) < FNU+F I M J /OF WCIM 
GAM3=-GAMJ*GAM2 
YNU = GAM2=«=BJ( 1 2)+YNU 
IF( AB.GT.O. )GO TO 120 
El=tXNU=^GAM2 
BY< 2)=E1^BJ( I2)+BY(2) 

IF( 12pGEpK)Gn TD 130 

E1^(GAM2-GAM3)/2p 

BYJ 2)=E1=«'BJ( 12+1 )+BY(2) 

120 GAM1=GAM2 
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133 
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135 
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140 

141 

142 
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121 GAM2=GAM3 140 

130 BYt 1 J =YNMJ 149 

IF( ABpGT.O. >BY< 2)={YNU*BJ(2)-TPX)/BJI 1 ) 150 


150 JAN=0 


151 


ABY = ABS< BY( 2) ) 

MAXs^NMl 

IF{ lNpE0,0)MAXt=K 
DO 160 I^ItMAX 
IF< JANpEOp I )GU T(j 160 
IMX==l+2 

1F( ABY.GTp U E6 1 JAN=1 

BY( 1+2 )==EAXc( J+FNU)'^BY( I + l )-BY< I ) 

ABY«ABS( BY( 1+2 ) ) 

160 CONTINUE 

IF( INpEO^O^GO to 300 
IP( IMXpLTpNPDGO TO 1 
BJJ( J)=BJ(NP1 ) 

BYY( J)=BY(NP1) 

1F{ JpGT. 1 IGO TO 210 
BJP(U=-BJ<2) 

BYP( U:=-By( 2) 

GO TO 220 
210 FAC=FNUP/X 

BJP( J )=-FAC*BJ(NPl ) + BJ(N) 

BYP( J)=-FAC<'BY( NPl )+BY( N) 

220 1F( JpLTp IOL )G0 TO 10 
N^J 

RETURN 

300 BJJ< l)=BJ( 1 ) 

BYY< U=:f3Y( 1 ) 

BJP( 1 )=-BJ( 2 ) 

BYP( 1 )=r-BY( 2 ) 

N = K 

IF( IMX.LTpK)N=IMX 
N=l+<N-l ) /INT 
1F«NpG7p1DL)N=I01 
DO 350 I=2|N 
L = l + ( I-1)=MNT 
LM=L-1 
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Fig. 8. 
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8JJn) = BJ(L> • 1®7 

8VY(I>=By<L) 188 

FAC=LM/X 189 

BJP( n=-FAC«BJ(L»+BJ(LM) 190 

350 BYP! I )=-FAC«BY( L >+BV{ LM) 191 

RETURN 192 

END 193 


Fig. 8. 


25 



SUBftaUTlwe SOROKC, S, I WR, 1 12f NEO) 
complex C( U^S( UtSS 

2 format ( IX, 115, lFl0*3WF15.7,lF10,0f 2F15.6) 

3 FORMATilHO) 

M = NEQ 

IFt Il2.ECl.2)r.G TQ 20 
C( 1 )=^CSQRT{Ct 1 ) ) 

00 A K=2,N 

A C(K)=Ci K)/C( 1 ) 

00 10 1 = 29 N 

lPn=!+l 

ID=( I-l 1*1-1 )/2 
II==IO+I 
00 5 L=l, IMG 
LI = K-1 )*N-{L*L-L )/2 + I 
5 C< I I)-C( I n-C( LI L I ) 

C( 1 1 )=CS0RTIC{ II)) 

IF( IPO*GT.N)GO TO 10 

00 8 J=IPO,N 

IJ=IO+J 

00 6 M=l, IMG 

M0=(M-1 )*N-t M*M-M)/2 

Ml=MD+l 

MJ=MO+J 

8 C( IJ)=C( IJ)-C(MJ)#CIMI ) 

8 C( IJ)=:C< IJ)/C( II ) 

10 CONTINUE 
20 S( 1)=S( D/U I ) 

DO 30 I=2rN 
IM0=I«1 

00 25 L = U1M0 

LI-( L-l )*N-(L*L~U/2 + I 

25 sn)«Si n--c(Li)*s(L) 

1 I = ( I-l )*N-( I*I-I )/2+I 
30 Si I )=S( 1 )/C( I I ) 

NN=nN+n*N)/2 

S(N)=S(N)/C(NN) 

NM0=N-1 ' 

DO 40 I-ltNMO 

K=sN-I 

KP0=K+1 

KD=< K-l )*N-(K*K-K)/2 
DO 35 L=KPO,N 
KL-KO+L 

35 S(K)=S{K)-C1KL)*S(L) 

KK = KD-**K 

40 S(K)sS(K)/C(KK) 

IF( IWR*L€*0) GO TO 100 
WRITE! 6,3) 

*CN0R = ,0 
DO 50 I=UN 
$A = CAflS!Sm) 

50 IF (SA*GT.CNOR IGNORES A 
1F(CN0R.LE*0. )CN0R-1. 

DO 60 1 = 1, N 
SS=S( I ) 

SA=CABS(SS) 

SNOR=SA/CNOR 
PH = .0 

IF! SA.GT.O* )PH = 57* 29578* AT AN2U1MAG( SS) tREAUSS) ) 
60 WR1TE(6,2) 1 f SNOR , 5A, PH , S S 
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Fig. 9. Subroutine SQROT. 
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WRITe(6,3) 0063 

100 RETURN 006<f 

END 0065 


Fig. 9. 
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SUBROUTINE GL J< LL , JJ t PHBR, RHCR , SGC . FL J I 0001 
DIMENSION SGCdSOl 0002 
DATA PI/3. 14159/ 0003 
J=JJ-1 0004 
SGJ = t-l )#>»J 0005 
*-=LL-l 0006 
FLJ=.5 0007 
GNU=L«PI/PHBR 0008 
SC=SGC(LL1 0009 
GNUS=6NU>»GNU 0010 
TEST = A8S(GNU-J=!=PI/PHCR) 0011 
OEN=GNUS-IJ«-PI/PHCR K»2 0012 
IF( TEST.GT. .001 IFL J = SGJ>kGNUS*SC/[)EN 0013 
return 0014 
END 0015 


Fig. 10. Subroutine GLJ. 


SUBROUTINE GN J ( M , J J , PHCR t SNC . FN J ) 0001 
dimension SNC(150I 0002 
DATA Pl/3. 14159/ 0003 
J=JJ-1 0004 
SGJ = ( -1 l^-^J 0005 
N=M-l 0006 
NS=N4N 0007 
SC=SNC(M) 0008 
FNJ=.5 0009 
TEST=ABS(N-J*PI/PHCR) 0010 
0EN=NS-( J*PI/PHCR)-*2 0011 
IF( TEST.GT. .001 IFNJ^SGJ^NS^SC/DEN 0012 
RETURN 0013 
END 0014 


Fig. n. Subroutine GNJ. 
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APPENDIX III 
SUBROUTINE SQROT 


This subroutine considers the matrix equation ZI = V which represents 
a system of simultaneous linear equations. If the square matrix Z is 
symmetric, SQROT is useful for obtaining the solution I with V given. 

NEQ denotes the number of simultaneous equations and the size of the 
matrix Z. 

On entry to SQROT, S is the excitation column V. On exit, the 
solution I is stored in S. Let Z(I,J) denote the synmetric square matrix. 
On entry to SQROT, the upper-right triangular portion of Z(I,J) is stored 
by rows in C(K) with 

(37) K = (I - 1)*NEQ - (1*1 - I) / 2 + J. 

If 112 = 1 , SQROT will transform the symmetric matrix into the auxiliary 
matrix (implicit inverse), store the result in C(K) and use the auxiliary 
matrix to solve the simultaneous equations. If 112 = 2, this indicates 
that C(K) already contains the auxiliary matrix. 

The transformation from the symmetric matrix to the auxiliary 
matrix is prograrmed above statement 10, and the solution of the simul- 
taneous equations is programmed in statements 20 to 40. If IWR is 
positive, the program below statement 40 will write the solution. 

SQROT uses the square root method described in Reference [4]. The 
original symmetric matrix Z and the upper triangular auxiliary matrix A 

are related by 

(38) Z = A’ A 

where A' is the transpose of A. 

The determinant of the symmetric matrix Z may be obtained by 
squaring the product of the diagonal elements in the auxiliary matrix. 

SQROT was developed by Dr. Robert G. Wickliff Jr,, now with Hewlett 
Packard, Colorado Springs, Colorado 80907. 
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